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Abstract 

Using the C* algebraic scattering approacli to study quasifree fermionic sys- 
tems out of equilibrium in quantum statistical mechanics, we construct the nonequi- 
librium steady state in the isotropic XY chain whose translation invariance has been 
broken by a local magnetization and analyze the asymptotic behavior of the expec- 
tation value for a class of spatial correlation observables in this state. The effect of 
the breaking of translation invariance is twofold. Mathematically, the finite rank per- 
turbation not only regularizes the scalar symbol of the invertible Toeplitz operator 
generating the leading order exponential decay but also gives rise to an additional 
trace class Hankel operator in the correlation determinant. Physically, in its decay 
rate, the nonequilibrium steady state exhibits a left mover-right mover structure 
affected by the scattering at the impurity. 

Mathematics Subject Classifications (2000) 46L60, 47B35, 82C10, 82C23 



1 Introduction 

In the mathematical study of open quantum systems, the role played by quasifree 
fermionic systems is an important one. Within the framework of algebraic quantum 
statistical mechanics, they not only allow for a powerful description by means of scat- 
tering theory on the one-particle Hilbert space over which the fermionic algebra of 
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observables is built, being thus ideally suited for rigorous analysis on many levels, but 
they also represent a class of systems which are indeed realized in nature, see, for 
example, Culvahouse et al. [16], D'lorio et al. [17], and Sologubenko et al. [24]. A 
special instance of this class is the finite XY spin chain introduced by Lieb et al. [20] 
and extended to the infinite two-sided discrete line by Araki [5] in the framework of 
C* -dynamical systems. As a matter of fact, this spin model can be mapped, in some 
precise sense, onto a gas of free fermions with the help of the Araki-Jordan-Wigner 
transformation. In order to study the effect of the breaking of translation invariance 
in this system, we choose the physically interesting and computationally convenient 
emptiness formation correlation observable. The so-called emptiness formation prob- 
ability (EFP), i.e. the expectation value of this observable in a given state, describes, 
in the spin picture, the probability that all spins in a string of a given length point down- 
wards. However, we would like to underline that the analysis is not limited to this 
observable but can rather be carried out for a broad class of spatial correlations. 

The asymptotic behavior of the EFP in the XY chain for large string length has 
already been analyzed for the cases where the state is a ground state or a thermal 
equilibrium state at positive temperature. In both cases, the EFP can be written as the 
determinant of the section of a Toeplitz operator with scalar symbol. Since the higher 
order asymptotics of a Toeplitz determinant is highly sensitive to the regularity of the 
symbol of the Toeplitz operator, the asymptotic behavior of the ground state EFP is 
qualitatively different in the so-called critical and noncritical regimes corresponding to 
certain values of the anisotropy and the exterior magnetic field of the XY chain, i.e., 
in ([T9]) below, the parameters 7 and A, respectively. It has been found that the EFP 
decays like a Gaussian in one of the critical regimes (with some additional explicit 
numerical prefactor and some power law prefactor), see Shiroishi etal. [23] and refer- 
ences therein. In a second critical regime and in all noncritical regimes, the EFP decays 
exponentially (in contrast to the noncritical regimes, there is an additional power law 
prefactor in the second critical regime whose exponent differs from the one in the first 
critical regime), see Abanov and Franchini [1 , 18]. These results have been derived by 
using powerful theorems of Szego, Widom, and Fisher-Hartwig, and the yet unproven 
Basor-Tracy conjecture and some of its extensions, see Widom [26] and Bottcher and 
Silbermann [Til [TSj- Furthermore, in thermal equilibrium at positive temperature, the 
EFP can again be shown to decay exponentially by using a theorem of Szego, see, for 
example, Shiroishi etal. [23] and Franchini and Abanov [18]. 

In contrast, out of equilibrium, the situation is more subtle. The typical open system 
consists of a confined sample which is coupled to extended ideal reservoirs at different 
temperatures. Using this paradigm, a translation invariant nonequilibrium steady state 
(NESS) has been constructed in Aschbacher and Pillet [11] for the XY chain using the 
scattering approach to algebraic quantum statistical mechanics developed by Ruelle 
(22] (for 7 = A = 0, this NESS has also been found by Araki and Ho [6] using a different 
method; moreover, using the latter approach, the magnetization profile at intermedi- 
ate but large times has been studied by Ogata [21]). In this NESS, the EFP can still 
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be recast into the form of a Toeplitz determinant, but now, tlie symbol is, in general, 
no longer scalar and regular. Due to the lack of control of higher order determinant 
asymptotics In Toeplitz theory with nontrlvlal Irregular block symbols, we started off by 
studying bounds on the leading asymptotic order for a class of general block Toeplitz 
determinants in Aschbacher (TJ. There, it turned out that suitable basic spectral in- 
formation on the density of the state is sufficient to derive a bound on the rate of the 
exponential decay of the EFP In general translation invariant fermionic quasifree states. 
This bound proved to be exact not only for the decay rates of the ground states and the 
equilibrium states at positive temperature treated in Abanov and Franchini [1 , 18] and 
Shiroishi etal. [23] but also for the translation invariant NESS In the isotropic XY chain 
analyzed in Aschbacher [9]. 

In the present paper, new results are obtained for the asymptotic behavior of a 
class of spatial correlations, and In particular, for the asymptotic behavior of the EFP, 
in a NESS of the isotropic XY chain whose translation invariance has been broken 
by a local magnetization In the form of a finite rank perturbation. Although such a 
spatial correlation can again be transformed from its initial Pafflan form into a scalar 
Toeplitz determinant, the effect of the breaking of translation invariance manifests Itself 
in a regularizatlon of the Toeplitz symbol and the appearance of an additional Hankel 
operator whose symbol is smooth. Hence, due to Peller's theorem, this operator Is 
of trace class, and the spatial asymptotics is governed by an exponential decay due 
to the invertlbility of the Toeplitz operator at nonvanlshing temperature. Moreover, the 
decay rate, determined by the Toeplitz symbol, exhibits the underlying left mover-right 
mover structure affected by the scattering at the Impurity (see also Aschbacher [8] and 
Aschbacher and Barbaroux [10J for the left mover-right mover structure of the NESS 
expectation of several other types of correlation observables). 

The paper is organized as follows. In Section [2| we set the stage for the nonequi- 
librium XY chain with impurity, construct Its NESS, and derive the basic expression for 
the NESS EFP Section |^then contains the asymptotic analysis of the NESS EFP Sev- 
eral ingredients of the proofs have been transferred to the appendix, as, for example, 
the construction of the wave operators by means of stationary scattering theory or the 
summary of the spectral properties of the so-called magnetic Hamiltonian. 



2 Nonequilibrium setting 

In this section, we will shortly summarize the setting for the system out of equilib- 
rium used In Aschbacher and Plllet [11]. In contradistinction to the presentation there, 
we skip the formulation of the two-sided XY chain as a spin system and rather focus 
directly on the underlying C*-dynamical system structure in terms of Bogollubov au- 
tomorphisms on a selfdual CAR algebra as In ArakI [5]. A C* -dynamical system Is a 
pair (21, r), where 21 Is a C* algebra and M 9 t r* g Aut(2t) a strongly continuous 
group of ^-automorphism of 21. For more Information on the algebraic approach to open 
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quantum systems, see, for example, Aschbacher etal. |T2]. 

For some given g N u {0}, the nonequilibrium configuration is set up by cutting 
the finite piece 

Zs ■= {x eZ \ - N < X < N} (1) 

out of the two-sided discrete line Z. This piece will play the role of the confined sample 
whereas the remaining parts, 

Zl := {x eZ\x <-{N + 1)}, (2) 
Zr := {xeZ\x>N + l}, (3) 

will act as infinitely extended thermal reservoirs, eventually carrying different tempera- 
tures, see Figure [1] 

The observables of the system are specified by the following selfdual CAR algebra 
over the wave functions on the chain. 

Definition 1 (Observables) Let 5^(f)) denote the fermionic Fock space over the one- 
particle Hilbert space of wave functions on the discrete line, 

\):=e(Z). (4) 

With the help of the creation and annihilation operators a*{f), a{f) e C{^{i))) with f e 
P) (where C{H) denotes the bounded linear operators on the Hilbert space H), the 
complex linear mapping i? : f)®^ ^ '^(S'(f))) 'S defined, forF ■= [/i, /2] e f)®^, by 

B{F):=a\h) + a{h). (5) 

The observables are described by the selfdual CAR algebra over f)®^ with antiunitary 
involution J generated by the operators B{F) e C{^{i))) for all F e f)®^, i.e. we have, 
forallF,Ge f)®l 

{B*{F),B{G)} = {F,G), (6) 
B*{F) = B{JF), (7) 

where JF -.= [/a, A] for all F -.= [/ija] e f)®^, the anticommutator of A, B g C{n) is 
{A, B} := AB + BA, and the scalar product in P)®^ is written as the one in f). We denote 
this algebra by'il -.= 2t([)®2, J). 

Remark 2 The concept of selfdual CAR algebras has been introduced and developed 
in Arak\ [3, 4J. Here, it is just a convenient way of worl<ing with the linear combination 
Also in view of future generalizations of the present paper, for example to the case 
of the truly anisotropic XY chain and other classes of correlations, we will sticl^ to this 
notation in the present context. 
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Figure 1 : The nonequilibrium setting for the XY chain. 



We next specify the Bogoliubov *-automorphisms on the selfdual CAR algebra 
which describe the time evolutions used for the construction of the NESS. 

Definition 3 (Dynamics) Let the coupling strength be k, > 0, and letu e £(()) be the 
translation given by {uf){x) := f{x - 1) for all f e i) and all x e Z. The XY, the 
decoupled, and the magnetic one-particle Hamiltonians h, ho, e C{^), respectively, 
are defined by 



h 
ho 



Re{u), 

h- {vl + vr), 

h + KV, 



(8) 

(9) 
(10) 



where the decoupling operators vl, vr e (with C^{H) the finite rank operators on 
H) and the operator v e C^{^) which breaks translation invariance have the form 



Vr 

V 



u 



-(N+l) 



Re( 



= Po- 



ol) 
(12) 
(13) 



Here, the projection pq e is given by po := {So, ■ )So, where Sr, e f) forx e Z denotes 
the Kronecker function (moreover, the real part of A e C{V.) is given by Rc{A) ■= 
{A + A*)/2). For all t e R, the XY, the decoupled, and the magnetic time evolutions 
are the Bogoliubov * -automorphisms t\tq,t^ e Aut(2l) defined on the generators 
B{F) e 21 with F e f)®^ by 



t\B{F)) 
r'oiBiF)) 



B{e''"°F), 
B{e''"^F), 



(14) 
(15) 
(16) 



where we setH -.^h® -h, Ho := ho © -ho, and Hb := /is © -hB- 



Remark 4 For the sake of an easy exposition, we restrict the analysis to the case k> 0, 
the case k <0 being strictly analogous. 

Remark 5 The magnetic Hamiltonian Hb e breaks translation invariance in the 

sense that the commutator [Hb,u ® u] = [hB,u] © -[hB,u] is nonvanlshing (where 
[A, B] := AB - BA is the commutator oi A,B e C{n)), i.e., for all / e I), it holds 



[/iB,H]/ = /^(/(-l)5o-/(0)(5i). 



(17) 
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Remark 6 Since Hb e anticommutes with the antiunitary involution J, the mag- 

netic Hamiltonian Hb generates a Bogoliubov transformation in the sense of Aral<i 
131 |4J, i.e. that, for all t g M, we have 



J] = 0. 



(18) 



The same also holds for the XY and the decoupled Hamiltonian H, Hq g £(f)^ 



Remark 7 As mentioned at the beginning of this section, this model has its origin in 
the XY spin chain whose formal Hamiltonian is given by 



^ = -i E {(1 + ^) -i^^-^'^ + (1 - ^) + } 



(19) 



where 7 g (-1, 1) denotes the anisotropy, A g M the external magnetic field, and the 
Pauli basis of C^""^ reads 



' 1 


" 




' 


1 " 




' -i " 




" 1 " 





1 


, 0-1 = 


1 







i 


, 0-3 = 


-1 



(20) 



The Hamiltonian h from ^ corresponds to the case of the isotropic XY chain without 
external magnetic field, i.e. to the case where 7 = and A = 0. 

The left and right reservoirs carry the inverse temperatures /3l and Pr, respectively. 
Pour fixer les idees, we assume, w.l.o.g., that they satisfy 



< /3l < < 00. 



(21) 



IVIoreover, for later use, we set (3 -.= + /3l)/2 and S := {(3r - (3l)/2. 

We next specify the state in which the system is prepared initially. It consists of 
a KIVIS state at the corresponding temperature for each reservoir, and, w.l.o.g., of the 
chaotic state for the sample. For the definition of fermionic quasifree states, see Ap- 
pendix |A| 

Definition 8 (Initial state) The initial state uq g is the quasifree state specified 
by the density So e of the form 



So := So - © So,+, 

where the operators so,± e £(f)) are defined by 

So,± := {I + e^'T\ 

and /co G £(f) ~ f)L © f)5 © is given by 



(22) 
(23) 
(24) 



Here, for a = L,S,R, we used the definitions := and := i*^hia g C{[)a), 

where io, : f)^ f) is the natural injection defined, for any f e f)„, byiaf{x) := f{x) if 
X eZa, and zero otheoA/ise. 
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Remark 9 Note that 5*0 g is well-defined, and that it satisfies the properties of a 

density given in Definition [32] of Appendix|A} 



Remark 10 The one-particle Hilbert space f) over Z = u u Z/j decomposes as 
f) ~ f)i © f}5 © It follows from ^ in Definition [3] that, w.r.t. this decomposition, 
the decoupled Hamiltonian ho does not couple the different subsystems to each other, 
indeed, i.e. we have ho = hi ® hs ® h^. 

As discussed in the Introduction, we pick the EFP correlation observable in order to 
study the effect of the breal^ing of translation invariance on nonequilibrium expectation 
values. This observable is defined as follows. 

Definition 11 (EFP) Letxo e Z andn e N. The EFP observable An g 21 is defined by 

2n 

An:=l[B{F,), (25) 

i=l 

where, for alii e N, the form factors g f)®^ are given by 

F2i-i := u'®u'Gi, (26) 
F2i := u'®u'G2, (27) 

and the initial form factors Gi, G2 g f)®^ look like 

Gi:= JG2:=[0,5.„-i]. (28) 

Moreover, the expectation value P : N [0, 1] of the EFP observable A„ g 21 in the 



NESS ujb g £{Ql) constructed in Theorem 17 below is denoted by 



P(n) ■.= coBiA„). (29) 



Remark 12 As for the name EFP, note that An = Hxlto and that, with Bn := 

nililto"^ a^., we have, for any state u e £:(2t), 

xo+n— 1 

< u{An) = u{B^B:) < \\Bnf < n P-f = l- (30) 

x=xo 

Remark 13 The analysis of this paper can also be carried out for different form factors. 
If we choose the initial form factors Gi =: [gi,i,gi,2\ g P)®^ for i = 1,2 to be of the 
completely localized form = auS^^^ for an g C and xu g Z with / = 1, 2, we cover the 
case Gi = [-6-1,6-1] and G2 = [60,60]. This choice describes the prominent spin-spin 
correlations afVf ^ see, for example, Aschbacher and Barbaroux [10J. 
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The following definition from Ruelle [22] introduces the concept of nonequilibrium 
steady state (NESS) in the framework of C* -dynamical systems. For the situation 
at hand, the C*-dynamical system is given in terms of the magnetic Bogoliubov *- 
automorphism group tb on the selfdual CAR algebra 21. 

Definition 14 (NESS) A NESS associated with the C* -dynamical system (21, tb) and 
the initial state uq e ^(21) is a weak-* limit point for T ^ oo of the net 



T > 



(31) 



Next, we define the time dependent correlation matrix of the EFP observable A„ e 21 
w.r.t. the initial state wo e £^(2t) and the magnetic dynamics g Aut(2t). 

Definition 15 (Correlation matrix) Let Fi g f)®^ /bri g N be the form factors of Defi- 
nition\Tl\ For all t g R, the skew-symmetric correlation matrix Vtn{t) g c^"^^" := {Ae 



A' 

its entries 



-A} (where A^ is the transpose of A) is defined, for all i,j 

-''^^JFi)B{e''^^Fj)), if i < j, 



, 2n, by 




ifi=j, 
if I > J. 



(32) 



Moreover, the matrices n^^{t),n^p{t),n^''{t),n^p{t) g c^^^^^n Qre defined, forij 
1, ...,2n andi < j, by 



u;o(5*(e''^nac(//B)JFi)5(e'*^Blac(//B)F,)), 
a;o(i?*(e^*''nac(^B)Ji^.)5(e'*^npp(//B)i^,)), 
u;o(5*(e^*^npp(i7B)Ji^05(e"^Blae(i^B)i^,)), 
u;o(5*(e^*^npp(/7B)JF^)5(e"^npp(i7B)F,)): 



(33) 
(34) 
(35) 
(36) 



and are to be completed as in (|32]) fori > j. Here, lac(-f^B), ipp(-f^B) e £(f)®^) are the 
spectral projections onto the absolutely continuous and the pure point subspaces of 
Hb, respectively 



The contributions which will play a role in the large time limit are defined as follows. 



Definition 16 (Asymptotic correlation matrix) Let Fi e f)®^ for i e be the form 
factors of Definitional^ The matrices Vl^, g c^"^^" are defined, for ij = 1, 2n 
andi < j, by 

:= iOo{B*{W{Ho,HB)JF,)B{W{Ho,HB)F,)), (37) 
fiPP := Yl oJoiB*iUHB)JF,)BiUHB)F,)), (38) 

eespec (Hb) 
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and are to be completed as in ([32) for i > j. Here, le(-f^B) e £(f)®^) denotes the 
spectral projection onto the eigenspace associated with the eigenvalue e in the set of 
eigenvalues speCpp(ifB) of H^, and the wave operator W{Ho, Hb) e is defined 

by 

W{Ho,Hb) := s - lime-^*^°e"^nac(i/B). (39) 

t—>-oo 

The following theorem establishes the existence and uniqueness of the NESS and 
yields an expression for the EFP in this NESS. 

From now on, whenever an entry of a sl^ew-sym metric matrix is written down, we 
always assume that the row index is strictly smaller than the column index. IVIoreover, 
for the definition of the Pfaffian, see Appendix |A} 

Theorem 17 (NESS and NESS EFP) There exists a unique quasifree NESS ub g 
Q(2t) associated with the C* -dynamical system (21, tb) and the initial state uq e £^(2t) 
whose density Sb g £([)®^) has the form 

Sb = W*{Ho,Hb)SoW{Ho,Hb)+ Y1 UHB)Sole{HB). (40) 

eespeCpp(_f/B) 

Moreover, the expectation value of the EFP observable in this NESS is given by 

P(r^)=pf(^]- + fiPP). (41) 

Proof We proceed similarly to the proof of Theorem 3.2 in Aschbacher et al. [13]. 
To this end, we note that the expectation value in the quasifree initial state uq e Q(2l) 
of the correlation observable A„ g 21 propagated in time with the magnetic dynamics 
Tb g Aut(2t) can be written, for all t g M, as the Pfaffian of the correlation matrix 



Quit) G C2"^2n f^Qp^ Definition 15 



u;o(r^(An))=pf(^^n(t)), (42) 

where we used ([T8]l in Remarl<|6]to commute the antiunitary involution J across the 
unitary group generated by Hb g In order to treat the argument of the Pfaffian, 

we mal<e use of assertion (a) in Theorem [36] of Appendix [B] which states that, for the 
singular continuous spectrum, we have 

spec,,(i/B) = 0. (43) 

Hence, injecting lac(-f^B) + Ipp(-f^B) = 1 g £([)®^) to the left of JFi and Fj in the corre- 
lation matrix entry 



"^«JFi,S'oe'*^«F,), (44) 
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the correlation matrix can be decomposed as 



(45) 



where the matrices on the r.h.s. of (|45]l are given in Definition [15} Since the NESS 
is constructed in the large time limit, we separately study this limit for all the terms in 
( |45l ). So, using that the initial state is invariant under the decoupled time evolution, i.e. 
[Ho, So] = 0, the first term can be written as 



-UHoJtHB^ 



(46) 



Thus, with the help of the Kato-Rosenblum theorem from scattering theory for pertur- 
bations of trace class type (see, for example, Yafaev [27]), we find 



lim fi^^(t) = n 



aa 
n 1 



(47) 



where we used that Ho - Hb e £°(f)®^), and the r.h.s. is given in Definition [T6} For the 
second term on the r.h.s. of (|45]l, we have the bound 



m{t)\ < \MHB)Soe^'''^lUHB)JFMFj\ 



(48) 



Since the pure point spectrum of Hb consists of the two simple eigenvalues iee, where 
cb is given in assertion (c) of Theorem [36] in Appendix [B} we have lpp(i/B) e 
and, hence, it follows that 



lim fi^P(t) = 0. 



(49) 



The same holds for riP^(t), of course. For the last term on the r.h.s. of ([45]l, using 

Ipp(^b) = Ic^Hb) + i^sAHb), we get 



J2 e-''^''-'\UHB)JFi, Sole'{HB)F,). 

e,e'G{±eB} 



(50) 



IVIoreover, since assertion (c) of Theorem [36] also states that ran (Ig^ {Hb)) c f) © and 
raii(l_eB(-ffB)) c © f), and since the density of the initial state 5*0 e £(f}®^) has the 
blocl< diagonal form given in ([22) of Definition [8| the terms in ([50) for different energies 
vanish, and, hence, the time dependence drops out of ([50). This leads to 



(51) 



for all t G M, where the r.h.s. is given in Definition [T6] Finally, since the Pfaffian 

pf : c^"x2" ^ c is a continuous mapping, we get 



lim — 

T^oo T 



T 



dt Wo(rB(^n)) 



lim — 



dt ^iiSlnit)) 



(52) 
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Note that we didn't make use of the specific structure of the form factors Fi. Hence, 
since the algebra of observables 21 is generated by the operators B{F) for F e f)®^, and 
since the mapping 21 9 A uo{t^{A)) g C is continuous uniformly in t g M, the relation 
( |52| ) defines the unique NESS ub g Q(2t). The form (|40]l of the density Sb follows from 
pT) ) and ([38]|. IVIoreover, due to the completeness of the wave operator and Remarl<[6| 
Sb has the defining properties of a density given in Definition |32] of AppendixjA} This is 
the assertion. □ 



3 NESS correlation asymptotics 



In order to approach the asymptotic behavior of the NESS EFP from Theorem 17 



we 



start off by studying more closely the two pieces of the asymptotic correlation matrix 



given in Definition [16} For this purpose, besides the position space, we will use the 
momentum space and the energy space defined before Theorem |35] of Appendix |A] 
and in Definition [39] of Appendix [C| respectively. 



Lemma 18 (ac-structure) The asymptotic correlation matrix il^^ g C 
composition 



2nx2n 
a 



has the de- 



(53) 



where the matrices QT'^ e c^"^^" 



)aa,+ 
^2i-12j 



are defined, for alii, j = l,...,n, by Qli:^^2j-i 
0, and the nonvanishing entries are given by 



Qaa,— 

^hi-12j 

Qaa,+ 



hB)Si+xg-i,S^W^{h, /Ib)5j+xo-i), 
{w+{h, /lB)5j+xo-l, S+W+{h, /lB)5i+xo-l)- 



(54) 
(55) 



Here, s± g are the density components of the translation invariant XY NESS given 
in Theorem\^ of Appendix^ and the wave operators w± {h, Kb) g are defined by 



w±{h, Hb) := s - lime"'^e-'*'^nac(/iE 



(56) 



where, from now on, all the spectral projections of Kb are denoted as the ones for Hb 
given in the Definitions\T^and\T^ with Hb replaced by Hb- 



Proof In order to rewrite the absolutely continuous contribution to the asymptotic 
correlation matrix from Definition [T6| we want to tal^e advantage of the fact that the 
operator 



S = W*{H^,H)SoW{Ho,H) 



(57) 
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is the known density of tlie translation invariant XY NESS {i.e. the NESS for k = 0) 
given in Theorem|35]of Appendix|A[ For this purpose, we use the chain rule W{Ho, Hb) = 
W{Ho,H)W{H,Hb) which is permissible since H - Ho,Hb- H e C^{H) (the wave 
operators W{Hq,H), W{H,Hb) e £(f)®^) are defined as in ([39]l with the appropriate 
replacements). Hence, the absolutely continuous contribution becomes 

= {W{H, Hb)JF„ SW{H, Hb)F,). (58) 

Using the block diagonal structure of the operators H, Hb, S e £([)®^) and plugging the 
explicit form of the form factors from Definition Qj] into ([58]l leads to the assertion. □ 



In order to evaluate the nonvanishing entries ([54]| and (|55]l from Lemma [T8 



we 

determine the action of the wave operators on completely localized wave functions. 
The main computations are carried out in Appendix |C} 

Proposition 19 (Wave operators) Letx e Z be any site. Then, in momentum space 
[) = L'^{T), the wave operators w±{h, Hb) e £([)) act on the completely localized wave 
function e f) as 

w±{h,hB)e^{k) = e^{k) .^'^li^^y. , (59) 

where we sete^{k) -.= S^{k) = e'^^ forallx e Z and forall k e 



Proof. Plugging | |141| | into ( |131| ) in Appendix [c] and applying f : ^ ^- ^ from ( |124| | to e^;, 
the action of the wave operator is expressed in energy space f) as 



w±ih,hB)Ue) = (27^)-l/2(l_e2)-l/4([(e + i^/^^)^(e-iv^ 



Vl — ± IK 



Applying f : P) -> ^ from ( |126| ) in Appendix[c]to ([60]l yields the assertion. □ 



Remark 20 The action ([59} relates to the action of the wave operator for the one- 
center ^-interaction on the continuous line by replacing sin(|A;|) by \k\, see, for example, 
Albeverio etal. ID. 

We next turn to the pure point contribution. 

Lemma 21 (pp-structure) The asymptotic correlation matrix Q^^ e c^"""^" has the de- 
composition 

fiPP = 5^1]PP^ (61) 
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where the matrices Q^^^^ e c^"^^^" ^re defined, for all i,j = l,...,n, by Qlf^^^^j-i '■= 
Qlf^ := ^lf2'j-i '■= ^2i-\2j -=0' ^'^Q' the nonvanishing entries are given by 

^lf-12j ■= (lpp(^B)'5i+xo-l,So,-lpp(/iB)5i+xo-l), (62) 
^2i2j-l '■= (lpp(^B)'^i+xo-l; ■5o,+ lpp(^B)'^i+xo-l)- (63) 

Here, so,± g £(f)) are the density components of the initial state given in Definition\^ 

Proof Using the block diagonal structures of H^, Sq e and plugging the explicit 

form of the form factors from Definition [lT| into (|38]) leads to the assertion. □ 



In order to evaluate the nonvanishing entries ([62) and (|63]l from Lemma [21] we 



determine the form of the projections onto the pure point subspaces of the magnetic 
Hamiltonian. A summary of its spectral properties is given in Appendix [B| 

Lemma 22 (Pure point projection) The projection onto the pure point subspace of 
the magnetic Hamiltonian satisfies 

dim(ran(lpp(/iB))) = 1, (64) 

and its range is spanned by an exponentially localized eigenfunction /b g f) of /ib e 
£(f)) with eigenvalue cb > 1. 

Proof See Theorem |36| in Appendix p[ □ 
Collecting the properties of the absolutely continuous and the pure point contribu- 



tions to the asymptotic correlation matrix from Lemma [T8] to Lemma [22} we get the 
following structural assertion. 

For the ingredients from Toeplitz theory referred to in the remainder of the present 
section, see, for example, Bottcher and Silbermann [14.15J. IVIoreover, we denote by 
Xa : -> {0, 1} the characteristic function of the set A c M. 

Proposition 23 (Determinantal structure) The NESS EFP is the determinant of the 
finite section of a Toeplitz operator, a Hankel operator, and an operator of finite rank. 
The symbol a e L°°(T) of the Toeplitz operator reads 

a = ipBS-,L + (1 - '^b)s~,r, (65) 
where the functions v?b, s±,„ e L°°(T) with a = L,R, are defined, for k e {-it, it], by 



s±Ak) ■■= -(l±tanh[i/3,cos(A;)]), (66) 

sin^fA;) 

see Figure\^ Moreover, the symbol of the Hankel operator is smooth. 



14 



W. H. Aschbacher 





Figure 2: The symbol a{k) > with k e (-vr, vr] for = 2, /J^ = |, and k = to the left 
and ft = I to the right. The nonvanishing magnetic field regularizes the symbol. 

Remark 24 In the limit k 0, we recover the symbol derived in Aschbacher [Z,^ for 
the translation invariant case. 



Remark 25 Note that for nonvanishing coupling, the characteristic function in (|67]l is 
smoothed out. This will play an essential role in the asymptotic analysis of the corre- 
sponding Toeplitz determinant. 



Proof. The total asymptotic correlation matrix, defined by Qn ■= + ^ C 



2nx2n 



has the 2 X 2 block substructure = [A 
defined, for i,j = 1, ...,n, by 



=1 ... „, where the matrices Aj, e C^""^ are 



h 





-hi 



if I < J, 

if i = j, 
if I > J, 



(68) 



and the entries are given by 



Qaa,— 
^hi-12j 
Qaa,+ 
^hi2j-l 



OPP,+ 



2j^ if « < J> 
if I < J. 



In order to rewrite the argument of the Pfaffian in P(n) = pf(ri„) from Theorem 17 



(69) 
(70) 

in a 



form more suited for the subsequent analysis, we want to apply a similarity transforma- 
tion to ri„. To this end, for any i,j = 1, 2n with i < j, we denote by i?'*^' g 0{2n) the 
elementary matrix whose left multiplication with any matrix A e c^"^^" exchanges the 
zth and jth row of A (where 0{n) stands for the orthogonal matrices in M"^"). Then, 
using the matrix R e 0(2rz) defined by R -.= /?[2(n-fc)+i,2{n-fc)+/+i]^ can 
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transform fi„ into off-diagonal block form, 



e„ 
-e* 



(71) 



where the matrix 0„ g C"""", called the reduced correlation matrix, Is defined by Its 
entries 6^^ := %, and, for all i,j e N, the numbers % are given by 



-Cji, if i > j. 



(72) 



Hence, using assertions (a) and (b) from Lemma|34]of Appendix |A} we get 

P(n) = pf(^]„) 



nfn— 1 ) 

(-l)^pf 
det(e„). 



e„ 
-0* 



(73) 



Let us next analyze the structure of 0„. In order to do so, we subdivide the discus- 
sion Into the following two cases w.r.t. the starting site xo e Z of the EFP string. 

Case 1: xo>0 

With the help of Lemma [44| of Appendix [D} we make the decomposition 9„ = Gt,™ + 



QH,n, where the matrix 9^,™ e C"""" has the entries QT,ij ■= dT,ij given by 6T,ij ■= bx^ij if 
i < j, and 6T,ij := -crji If i > j- Here, for all i,j e N, we define 



bT,ij 

-CT,ji 



{ei-j, s-eo) + (e. 



a_eo), \U<j, 
[ej-i, s+eo) + {ei-j, a+eo), if i > j, 



(74) 
(75) 



where s±,a± e L°°(T) are given In Theorem [35] of Appendix |A] and Lemma [44] of 
Appendix [□} respectively. Similarly, the matrix O^,™ e C"""" has the entries QH,ij ■= 
eH,ij given by 0^,^ := if i < j, and ^^.y := -CHjt if i > j- Here, for all i,j e N, we 
define 



-CHJi 



{ei+j,b_eo), \U<j, 
(ei+j,6+eo), if«>j, 



(76) 
(77) 



where 6± g L°^(T) are given in Lemma|44]of Appendix [D| Using ([66]l and ([67]| in (74^ , 
and - = - s„,l) and (ei_j, s+eo) = -(e^-i, s_eo) in ([75), where the last 

identity Is due to Definition [32] of the density In Appendix [A] we have bT,ij = (ei_j, aeo) 



for i < j and -crji = {ei-j, aeo) for i> j. Hence, O^.n is the finite section of the Toeplitz 
operator T[a] e C{f{N)) generated by the symbol a e L°°(T), i.e. 



QT,n = Tn[a]. 



(78) 
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Moreover, as for (|76]) and (|77]), using ( |1 19[ ) in Remark |37] of Appendix [B} we find that 

bH,ij = (ei+j_i,6eo) for i < j and -cnji = (ej+j_i,6eo) for i > j, where the function 
b e C°°(T) is defined by 



b{k) := IK, 



-ik{2xo-l) 



sin(fc) + i/t [ 



{fB,So-fB) 



(79) 



Hence, Qn-a is the finite section of the Hanl<el operator H[b] e C{f(N)) generated by 
the symbol 6 e C^iT), i.e. 



QH,n = Hn[b]. 



(80) 



Therefore, it follows from (|73]) that the NESS EFP is the determinant of the finite section 
of the sum of a Toeplitz and a Hanl^el operator. 



P(n) = det(T„[a]+i/„[6]), 



(81) 



where, in this case, the finite rank operator from the formulation of the assertion van- 
ishes. 

Let us now turn to the case where the EFP string starts to the left of the origin. 
Case 2: xq < 



For n > 1 + no, where we set uq ■= -xq, we again have from Lemma 44, that, for all 



,n-no, 



On. 



i+noj+no 



(82) 



where we set c := e_(2„„+i)fe. Defining the operator : C"« © C"« © e^{N) on 

/] G C"" © f{N) by the matrix multiplication with the infinite matrix 6ij from ([72]|, we 
have 



M := e - © (T[a] + H[c]) G £°(C"° © 



(83) 



Since the reduced correlation matrix satisfies 6„ = -Rn0-Rntran(ii„) with i?„ := 1 © P„ g 
£(C"o © fiN)), it follows from ([73| that 

P(n) = det(0 © (T„_„Ja] + Hn-no[c]) + Mn), (84) 

where M„ := RnMRn\ri,n{R„)- Hence, we arrive at the assertion. □ 



We are now ready to formulate our main result on the behavior of the NESS EFP 
for large string lengths. 

Theorem 26 (Exponential decay) Forn oo, the NESS EFP P : N ^ [0, 1] has an 

exponentially decaying bound, 



F{n) = 0{e 



-rn\ 



(85) 
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The decay rate r := r^? + Tb > contains the two parts 

Tr ■■= -t; tt ^og[s-,Rm (86) 
1 n c\k 

where the function (Tb g L'^ (T) is given by 

sm (/e) + 



Remark 27 Note that Theorem |26| holds for any coupling strength. In the small cou- 
pling limit, we recover the exact decay rate from Aschbacher [9J, and the bound derived 
for general quasifree systems In Aschbacher [7J. 

Remark 28 As can be seen in ([87]|, the NESS EFP decay rate displays a left mover - 
right mover structure. It is composed of a left mover carrying temperature Pr and com- 
ing from +00, a left mover carrying temperature Pn having been reflected at the pertur- 
bation at the origin, and a right mover carrying temperature /3l having been transmitted 
through the origin. This left mover - right mover structure has already been observed In 
translation invariant systems for several types of correlation functions, see Aschbacher 
(HUD and Aschbacher and Barbaroux [10J. 



Remark 29 Defining analogously to ([86]l, we have the rewrltings, for a = L,R, 

/f At. 
— log[i(l-tanh2[^/3«cos(A;)])], (89) 
'2 

/f At. 
— log [1(1 - [(1 - aB{k))tanh[ll3ncos{k)] 
_K It: 



+(TB(fc)tanh[|/3iCos(fc)]]2)l. (90) 



From ([89) and ([90]l, It immediately follows that, if the system is truly out of equilibrium, 
i.e. if 5 > 0, the decay rates are ordered as 

o<rz.<rB<rH, (9i) 

see also Figure [3} 

Remark 30 It follows from assertion (a) In Proposition [45] of Appendix[E]that the nonva- 
nlshlng coupling regularizes the underlying Toeplltz theory In the sense that the symbol 
which determines the decay rate is smoother than in the case k = 0, see Figure [2} 
Namely, the latter case requires FIsher-HartwIg theory and. If 5 > 0, leads to a strictly 
positive power law subleading order as given In Aschbacher [9]. 
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Proof. Since the Toeplitz symbol a e L°°(T) from ([65]l is real-valued, we can make 
use of the Hartman-Wlntner theorem in order to control the spectrum of the selfadjoint 
Toeplitz operator T[a] e C{f{N)). Moreover, due to Proposition |45] of Appendix [E| we 
have 

a e C(T), (92) 

and, hence, spec(T[a]) = ran (a) = [s__/j,(0), s+,ij(0)], where < s_,/j(0) < s+,r(0) < 1 in 
the temperature range < (3l < Pr < oo. Therefore, T[a] e C{f{K)) is invertible, 

i spec(T[a]), (93) 

and the spectrum is independent of the coupling strength. IVIoreover, since ([92) and 
( |93| ) hold, the Gohberg-Feldman theorem implies that the sequence {T„[a]}„6N is sta- 
ble, 

limsup||r-^[a]|| < oo. (94) 

For the following analysis, as in the proof of Proposition [gSj we discuss the cases 
xo > and xo < separately. For convenience of exposition, we start with the second 
case. 

Case 2: xq < 

Since we want to analyze the asymptotic behavior for large n with the help of Szego's 
strong limit theorem, we write, using | |94| |, 

P(^) ^ P(n) det(T„_„Ja]) det (T„[a]) 

G{aY det(T„_„„[a]) det(T„[a]) ^(a)" ' ^ ^ 

where G{a) is the exponential of the 0th Fourier coefficient of log(a), and no := -xq as 
before. Due to ([84]l and (|94]), the first factor on the r.h.s. of ([95]l can be written as 

'Pin] 

^^^^rj}^ = det(l + 1 © T-i„Ja]((-l) © iJ„_„Jc] + M„)). (96) 
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Moreover, since we know from Proposition |23] that c e C°°(T), we also have 

cG L°^(T) n5|(T), (97) 

where 5^ (T) are the usual Besov spaces. Therefore, Peller's theorem allows us to 
conclude that [C^iH) are the trace class operators on the Hilbert space U) 

H[b] G C\e\N)). (98) 

Due to ([83]l, (|84]), and ([98]l, the r.h.s. of ([96]l converges to the constant K{a,b) := 
det(l + 1 © T-^[a]{-l © H[b] + M)). In order to treat the second factor on the r.h.s. 
of ([95]l, we apply Szego's first limit theorem which is applicable due to (|92]) and ( [93] ). 
Hence, if we factorize the quotient as 

det(T„_„Ja]) ^ det (T„_4a]) 



det (T„ [a] ) ^J- det (T„+i_i [a] ) 



each factor on the r.h.s. of ([99]l converges to l/G{a). In order to treat the third factor 
on the r.h.s. of (|95]), we mal^e use of Szego's strong limit theorem. This theorem states 
that, since a{t) > for all t e T s.t. ind(a) = 0, and since 

aeW{T)r\Bl^\T), (100) 



which follows from a e C^(T) n PC°°(T) in Proposition |45| of Appendix [E| (vr(T) is 



the Wiener algebra and PC°°{T) are the piecewise smooth functions), the quotient 
converges to a constant usually denoted by E{a). Plugging the foregoing three limits 
into the r.h.s. of ([95]l, we get 

P(n) 

\im ^^ = K{a,b)Eia)Giar\ (101) 



Hence, in order to determine an exponential bound on the asymptotic decay, we are left 
with the computation of the constant G{a). Decomposing the integral in the 0th Fourier 
coefficient of log(a) w.r.t. positive and negative momenta and using the fact that s^,^ 
for a = L,R and ctb are even functions in A; g {-it, tt], we arrive at the expressions 
and (|87]) for the decay rate of the bound on the exponential decay. 

The case where the EFP string starts at nonnegative sites is simpler and is treated 
analogously as follows. 

Case 1: xo > 

Writing ([81} as in ([95]l, where, in this case, the second factor is absent, we can proceed 
as for Case 2. In particular, the determinant of the Toeplitz contribution can again be 
separated due to ([94), ([98]l holds for the symbol b satisfying ([97]|, and the constant now 
reads K{a,c) = det(l + T~^[a]H[c\). Finally, the last factor in dlOlj ) is absent. Hence, 



we arrive at the assertion. □ 
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Remark 31 The study of the present problem for the anisotropic XY model, i.e. for the 
case where 7 7^ in ([T9]| of Remarl< |7| is more complicated. Not only the Pfaffian 
structure of the correlation cannot be preserved in the present form, but also one has 
to cope with Toeplitz theory for operators with nonscalar symbols. We will study this 
set of problems for general quasifree systems elsewhere. 



A Fermionic quasifree states 

Let 21 be the selfdual CAR algebra from Definition [I] We denote by £:(2t) the set of 
states, i.e. the normalized positive linear functionals on the C* algebra 21. 

Definition 32 (Density) Tine density of a state u e £^(2t) is defined to be tfie operator 

S e witti 0<S* = S<land JSJ = 1-5 satisfying, for all F, G e f)®^ 

uj{B*{F)B{G)) = {F,SG). (102) 
An important class are the quasifree states. 

Definition 33 (Quasifree state) A state u e £^(2t) is called quasifree if it vanishes on 
the odd polynomials in the generators and if it is a Pfaffian on the even polynomials in 
the generators, i.e. if, for all Fi,..., F2„ g f}®^ and for any n eN, we have 

a;(5(Fi)...5(F2„)) =pf(fi„), (103) 

where the skew-symmetric matrix Qn e C^''''^" = {A e c^"""^" | A* = -A} is defined, 
forij = l,...,2n, by 

(u{B{F,)B{F^)), ifi<3, 
%:=<^0, if 1 = 3, (104) 

[-u{B{F,)B{Fi)), ifi>j. 

Here, the Pfaffian pf : _^ ^ ,5 gj^Qf^ 

n 

Pf(^) Xl^^Sn(vr) JjA^(2j-i),^(2j), (105) 

TT j = l 

where the sum is running over all pairings of the set {1,2, ...,2n}, i.e. over all the 
{2n)\/{2'"-n\) permutations n in the permutation group of 2n elements which satisfy 
7r(2j - 1) < 7r(2j + 1) and 7r(2j - 1) < 7r(2j), see Figure^ The set of quasifree 
states is denoted by Q(2l). 



The following lemma has been used in Section |3} 
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nnn-nrTn-nrm-iTTin-- 

7T = (123456) 71 = (123546) vr = (123645) vr = (132456) 

Figure 4: Some of the pairings for n = 3. The total number of intersections / relates 
to the signature as sign(7r) = (-1)^. 

Lemma 34 (Pfaff ian) The Pfaffian has the following properties. 

(a) LetX,Y e c^^^^n withY' = -Y. Then, 

pf (xrx*) = det(x) pf (y). (1 06) 

(b) LetX e C"^". Then, 

)n(n-l) 

= (-1)^ det(X). (107) 

Proof. See, for example, Stembridge [251. n 

Next, we state the properties of the NESS for the translation invariant case k = 0, 
the so-called XY NESS. To this end, let f : f) ^ ^ := ^^(T) (with unit circle T) 
be the Fourier transformation defined with the sign convention f{k) -.= {ff){k) -.= 
^^gg/(x)e'''^. IVloreover, for any a e £(P)), we use the notation a := faf*. We then 
have the following. 

Theorem 35 (XY NESS) There exists a unique quasifree NESSu e Q(2t) associated 
with the C* -dynamical system (21, r) and the initial state uq e Q(2l) whose density 
S E /:(f)®2) has the form 

S = s_®s+, (108) 
where the operators s± e C{i)) act in momentum space ^ as multiplication by 

s^{k):=^-{l±g^{k)), (109) 

and the functions g± : T ^ [-1,1) are defined by 

g±{k) := tanh[i(/3±sign(sin(A;))5)cos(fc)]. (110) 



pf 



X 
-X* 



Proof See Aschbacher and Pillet [TT]. 



□ 
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B Magnetic Hamiltonian 

In this section, we summarize the spectral theory of Hb e £(f}®^) needed above. To 
this end, we denote by specsc(^), spec^^{A), and speCpp(A) the singular continuous, the 
absolutely continuous, and the pure point spectrum of the operator A, respectively. 

Theorem 36 (Magnetic spectrum) The magnetic Hamiltonian e £(P)®^) has the 

following properties. 

(a) spec,,(i/B) = 

(b) spec,,(i/B) = [-l,l] 

(a) speCpp(ifB) = {icB} witheB > 1 
The eigenvalues ±eB are simple, and 

eB = VlT^. (111) 

The normalized eigenfunction of Hb with eigenvalue cb is given foy /b © e f)®^, 
where /b is exponentially localized, i.e. for all x e Z, it has the form 

/b(x):=— e-^BN^ (112) 
riB 

and the decay rate and the normalization constant are given by 

Ab := log(fi: + eB), (113) 

(114) 

Moreover, the eigenfunction of Hb with eigenvalue -cb reads © /b g f)®^. 

Proof Assertions (a) and (b) are proven in a more general context in Hume and 
Robinson [19] (which also contains the case of the truly anisotropic XY model without 
magnetic field, and more general perturbations). The fact that there is no eigenvalue 
embedded in the continuum speCj,c(i7B) = spec^^(iJ) = spec(iJ) = [-1,1] also follows 
from [19]. Hence, in order to derive assertion (c), we compute eigenfunctions of the 
operator Hb by loo[<ing for solutions of the eigenvalue equation /ib/ = e/ for e g M with 
|e| > 1 and not identically vanishing f e\). Since such e lie in the resolvent set of h, we 
can write / = -{h - e^^vf = -Kf{0){h - e)"^5o- By taking the scalar product of this 
equation with 5^, for any x g Z, we have 

f(x) = -KfiO)i6^,,ih-e)-%), (115) 

implying that /(O) ^ 0. Plugging x = into ( |1 15| ), we get the eigenvalue equation 




1 + k{So, {h - e)-i(5o) = 0. 



(116) 
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Switching to tlie momentum space representation and using Caucliy's residue theo- 
rem, we get, for all e g M with |e| > 1 and all x eZ, 



(5,, {h - e)-'6o) 



si: 



-sign(ej 



ign(e)\/e^ 



1) 



\x\ 



(117) 



If we plug X = into ( |1 17| | and use the assumption k > 0, we see that the equation 
dl 16| ) can be satisfied for e > 1 only. Solving ( |1 16| ) for this case leads to ( |1 1 1| ). Next, 
plugging ( pTl) into ^JT}, we get /(x) = /(0)(/€+eB)"l'^l from ( pTst . Choosing /(O) > 0, 
we arrive at ( |112D with | |1 13| and ([T]^- □ 



Remark 37 The Fourier transformation of /b g f) is given, for all k e (-vr, tt], by 



fBik) 



n-B cb — cos(fc) 

Cauchy's residue theorem yields that, for a; g Z with a; > 0, we also have 



(118) 



dk 



-ikx 



^ 271 sin(A;) + m 



(119) 



The integrand on the r.h.s. of (119) is used for the extraction of the Hankel symbol in 



the proof of Proposition 23 



C Wave operators 

In this section, we use the stationary approach to scattering theory in order to compute 
the wave operators w±{h, Hb) g £(f}) appearing in the ac-contribution to the asymptotic 
correlation matrix from Lemma [18} To this end, we first need to express the resolvent 
of the magnetic Hamiltonian by the resolvent of the XY Hamiltonian. This is done in the 
following lemma. 

For any operator a g £(f)) and any 2 g C in the resolvent set of a, we denote by 

r^{a) := {a - zy^ G £(f)) the resolvent of a at the point z. 

Lemma 38 (Magnetic resolvent) Lete eR ande > 0. Then, at the points e ± ie, the 
resolvent of e C{[)) can be expressed in terms of the resolvent of h e C{[)) as 

re±\e{hB) = re±ie{h) - — — (re=pie(/l)5o, " ) re±ie{h)6o. (1 20) 

1 + K{do,re±ie{h)do) 

Proof In order to simplify notation, we drop the index e±ie oi the resolvents. With the 
help of the resolvent identity r(/iB) = r{h) - Kr{h)vr{hB), we can write, for all / g P), 



r{hB)f + Ki6o, r{hB)f)rih)6o = rih)f. 



(121) 
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Taking the scalar product of ( |121| | from the left with So, we get 

(l + «:(5o,r(/i)5o))(^o,r(/iB)/) = (r(/i)*5o, /)• (122) 

Since, due to (1 + K,{6o,r{h)6o)){l - K{6o,r{hB)6o)) = 1, the first factor on the l.h.s. 
of is nonvanishing, we can solve ( |122| for {6o,r{hB)f). Plugging the resulting 
expression into ( |121[ ) yields the assertion. □ 



In the next definition, we introduce the energy space being the direct integral de- 
composition of the absolutely continuous subspace of the XY Hamiltonian h e £([)) 
w.r.t. which h is diagonal. 

Definition 39 (Energy space) Let the direct integral over spec^^{h) with fiber be 
denoted by 

t):=L\[~l,l],C'), (123) 

and let us call ^ the energy space of h. Moreover, the mapping ~f -.[) ^ [) is defined, for 

all if e [), by 

{fip){e) := (27r)-^/2(l -e')-^/'^[v9(arccos(e)),(^(-arccos(e))]. (124) 

We will use the notation f -.= fff for all f Ei),anda = ffaf f for all a e C^)), where the 
Fourier transform f -.[)-)■[) is defined in Appendix^ Moreover, the Euclidean scalar 
product in will be denoted as {■,■). 

We then have the following lemma. 

Lemma 40 (Diagonalization) The mapping } e C{[), [)) is unitary and the XY Hamilto- 
nian h G £([)) acts, on any r] ei), as the multiplication by the energy variable e, 

Chr]){e) = er]{e). (125) 

Proof A simple computation shows that f is a surjective isometry with f ^ = f : ^ ^ 
acting on all 77 =: [r]i,r]2] g P) as 

ifr])ik) = (27r)i/2(i _ cos2(fc))V4 [x[o,^](A;) r7i(cos(fc)) + X[~.,o]ik) V2icosik))]. (126) 
Equality ( |125| then follows immediately. □ 

We introduce the following abbreviations. 
Definition 41 (Boundary values) Lete g M ande > 0. For all f,g g f), we define 

Q},g{e±ie) := [f .re±ie{h)g), (127) 
7/,3(e,^) := ^(^/,9(e + ie) -^/,5(e-ie)). (128) 
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Moreover, if the limits exist, we write 



^/,9(e±iO) := lim ^/,g(e±ie), 

£-5>0+ 

7/,5(e) := lim 7/,5(e,e). 

e-s>0+ 



(129) 
(130) 



The wave operators then have the following form. 

Proposition 42 (Wave operators) In energy space t), the action of the wave operators 

w±{h, /ib) g £(f)) on any / g f) has the form 



w^{hM)he) = /(e) - .l^'^'^'flL ^o(e). 

1 + f^QSoMe ± lO) 



(131) 



Proof In order to compute the wave operators w±{h,hB) e £(P)) with the help of the 
stationary scheme in scattering theory (see, for example, Yafaev [27J), we write them 
in the weak abelian form 



w±{h, hB) = w- lim 2e I dt e-''''l^{h)e^''''e'^'"'^U,{hi 



(132) 



Applying Parseval's identity to ( [132) and using that r^±i,{h) = ±i dt eTi*(''-(<=±i^)), we 
can write, for a\\ f,g ei). 



if,w±{h,hB)g) = lim - 

£-s>0+ TT 



de (r {h)f, re±ie{hB)lac{hB)g)- 



(133) 



Moreover, if the limits e 0+ of ire±ie{h)f,re±ie{hB)g) exist for all f,g e i) and almost 
all e G M (the set of full measure depending on / and g), we get 



(/, w±{h, hB)g) 



I de lim -{re±\e{h)f ,re±ie{hB)g) 

J_l £^0+ TT 



(134) 



because Is^dh) = 1 and spec{h) = [-1, 1]. In order to compute the limit in | |134| |, we 
express the resolvents re±ie{h) of the magnetic Hamiltonian in terms of the resolvents 
re±ie{h) of the XY Hamiltonian. Plugging from Lemma [SS] int o the scalar product 
on the r.h.s. of and using and ( |128[ ) from Definitiorr|4T] we have 



- {re±ie{h)f,re±ie{hB)g) = 7/,9(e,^) - 7— / , . X 

TT 1 + Kgso,SQ{e ± le) 



lf,SQ{e,e) gso,gie±ie). (135) 



Here, we made use of the fact that, due to the resolvent identity, we have the equality 

7/,g(e,£:) = (/, ^re±ie{h)rezpie{h)g). Now, we know that, for any f,g e [) and almost all 
e G [-1, 1], the following limits exist. 



^/,g(e±iO) = ± 



.dif,pie)g) 



■Ri- 



de 



\, 1 d{f,p{e')g) 



de' 



(136) 
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where the p.v.-integral denotes Cauchy's principle value, the mapping p B{R) £([)) 
with B{R) the Borel sets on R is the projection-valued spectral measure of the XY 
Hamiltonian h, and we used that 

d{f,p{e)9) 



IVIoreover, we get from ( |128D and ( |136D , 



de 



de. 



7/,s(e) 



d{f,p{e)g) 
de 



Therefore, plugging ( |135| ), ( |136| ), and ( |138[ ) into ( |134| |, we can write 

if,w±{h, h^)g) = {f,g) - k 



^ 7/.3o(e)g5o,g(e±i0) 
_i 1 + K^5o,5o(e ± iO) ' 



(137) 



(138) 



(139) 



where, in the first term on the r.h.s., we used X^^de 7/,g(e) = (/, li,c{h)g) = U',9)- '"^ 



order to write the derivatives in ( |138[ ) entering ( |139| ) more explicitly, we switch to the 



energy space representation from Definition |39 
have, for a\\ f,g ei), that 

d(f,p{e)g) 



Using the diagonalization ( |125| l, we 



de 



(/(e), ^(e)), 



(140) 



where we recall from Definition 39 that 
of the direct integral ^ 



) denotes the scalar product in the fiber 
L2([-i;i],C2), and / = fff for all f e [). Hence, plugging 
( |138| and ( |140| | into ( |139[ ), we arrive at the assertion. □ 



Finally, since the wave operators w±{h, h^) e £(f}) appearing in the ac-contribution 
to the asymptotic correlation matrix act on completely localized wave functions 4 g f) 



with X G Z, we compute the terms gso,5, (e ± iO) on the r.h.s. of ( |131[ ) in Proposition |42 

Lemma 43 (Boundary values) Letx g Z ande g (-1, 1). Then, we have 

(e^iv/r^)!^! 



^<5oA(e±iO) = ±i 



(141) 



Proof. Letx g Z with x > 0, e g (-1, 1), and e > sufficiently small. Writing gsa^s^e-ie) 
in the momentum space representation, using Cauchy's residue theorem, and taking 
the limit e -> 0+, we get the expression ( |141[ ) for gso^i^^ - iO)- IVIoreover, using ( |127| |, 
the translation and parity invariance of h, i.e. [h,u] = and [h,6] = 0, respectively, 
where 6 : i) ^ [) \s defined, for all / g P), by {9f){x) := /(-x), we have, for all x g Z 
with X > 0, that 



P5o,5.(e + ie) 
^<5o,5-.(e - ie) 



This yields the assertion. 



(142) 
(143) 

□ 
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D Asymptotic correlation matrix 

In this section, we compute the nonvanishing entries of the total asymptotic correlation 
matrix used above. 

Lemma 44 (Structure) Letxo > andij > 1, orxo < andij > 1 - xq. Then, the 
entries of the asymptotic correlation matrix have the structure 

^2i-i2j + ^2i-i2j = (ei-j, s-eo) + {et-j, a_eo) + (ei+j, 6_eo), ifi < j, (144) 

)aa,+ I QPP,+ 

where the functions a±,b± e L°^(T) are defined by 



^2?2t-i + ^2%-! = (ej-i, s+eo) + (ei_j, a+eo) + (6^+^, 6+eo), ifi > j, (145) 



a,(fc) := - Xio,.m—^^^y^^^ (146) 

S, r-rt"! *<■ „-2AB(a;()-l) 

b±{k) = (-i^) e-^'^(-°-^) + 4(/B,.o,^/B)- ^. (147) 

sm(A;) + i/t rig cb — cos(A;) 

Proof For z < j, plugging the wave operator (|59]) into ac-contribution (|54]| yields 

d)^ g-i(fc(i+a;o-l)-|fc||i+xo-l|) 



sin(|A;|) — in 



dk g-i(|fc||i+xo-l|-A:(i+a;o-l)) 

^^^'-^^'^ sin(|A:|) + i«: 

dk „-i\k\{\i+xo-l\-\j+xo-l\) 



Moreover, using Lemma 22 we have 



Q^ll,^ = ^(/B,so,-/B)e-^B(l^+^'°-^l+l^+^»-^l). (149) 

The expressions for i > j are analogous. Then, if xq > and i,j > 1, or xq < and 
i,j > i - xo, we use the translation invariance of s±, resolve the absolute values, and 
decompose the integrals w.r.t. the sign of the momentum in order to get rid of the sign 



function in the density s±. This leads to (144 1 and (145|. □ 



E Toeplitz symbol regularity 



The following proposition is used in Theorem 
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Proposition 45 (Regularity) The Toeplitz symbol a e L°°{T) of Proposition\2^ has the 
following properties. 



(a) a e Ci(T) nPC°°(T) 

(b) The left and right derivatives D±a'{k) exist for all k e (-vr, n], but, for 
k- := vr, we have 



D+a'{k±)-D^a{k±] 



1 smh[|(/3^ - /jz.)] 
cosh[^/3/j] cosh[|/3i] 



and 



(150) 



Proof. From the very form of the symbol given in (|65]) - (|67]|, we get a e PC°°{T) with 
jumps at k±, and since ips g C(T) for nonvanishing coupling, we also have a e C{T). 
IVIoreover, the one-sided limits yield a'(A;± + 0) = a'(A;± - o) = which is assertion (a), 
and analogously for assertion (b). □ 



Acltnowledgements We would lil<e to thanl< the editor and the referee for their con- 
structive remarl^s. IVIoreover, the support provided by the German Research Founda- 
tion (DFG) is gratefully acl^nowledged. 



References 

[1] Abanov G A and Franchini F 2003 Emptiness formation probability for the 
anisotropic XY spin chain in a magnetic field Phys. Lett. /\ 316 342-9 

[2] Albeverio S, Gesztesy F, Hoegh-Krohn R, and Holden H 2000 Solvable models in 
quantum mechanics (AlVIS: Providence) 

[3] Aral^i H 1968 On the diagonalization of a bilinear Hamiltonian by a Bogoliubov 
transformation Publ. RIMS Kyoto Univ. 6 385-442 

[4] Aral<i H 1971 On quasifree states of CAR and Bogoliubov automorphisms Publ. 
RIMS Kyoto Univ 6 385-442 

[5] Aral^i H 1984 On the XY-model on two-sided infinite chain Publ. RIMS Kyoto Univ 
20 277-96 

[6] Aral^i H and Ho T G 2000 Asymptotic time evolution of a partitioned infinite two- 
sided isotropic XY-chain Proc. Steklov Inst. Math. 228 191-204 

[7] Aschbacher W H 2007 On the emptiness formation probability in quasi-free states 
Cont. Math. 447 1-16 

[8] Aschbacher W H 2007 Non-zero entropy density in the XY chain out of equilibrium 
Lett. Math. Phys. 79 1-16 



Broken translation invariance in quasifree nonequilibrium correlations 



29 



[9] Aschbacher W H 2010 A remark on the subleading order in the asymptotics of 
the nonequilibrium emptiness formation probability Confluentes Math. 2 293-31 1 
(arXiv:1 009.1 584J 

[10] Aschbacher W H and Barbaroux J IVI 2007 Exponential spatial decay of spin-spin 
correlations in translation invariant quasifree states J. Math. Phys. 48 1 13302-1 - 
14 

[11] Aschbacher W H and Pillet C A 2003 Non-equilibrium steady states of the XY 
chain J. Stat. Phys. 112 1 153-75 

[1 2] Aschbacher W H, Jaksic V, Pautrat Y, and Pillet C A 2006 Topics in non-equilibrium 
quantum statistical mechanics Lecture Notes in Mathematics 1882 (New York: 
Springer) 1-66 

[13] Aschbacher W H, Jaksic V, Pautrat Y, and Pillet C A 2007 Transport properties of 
quasi-free fermions J. Math. Phys. 48 032101-1-28 

[14] Bottcher A and Silbermann B 1999 Introduction to large truncated Toeplitz matri- 
ces (New York: Springer) 

[15] Bottcher A and Silbermann B 2006 Analysis of Toeplitz operators (Berlin: 
Springer) 

[16] Culvahouse J W, Schinke D P, and Pfortmiller L G 1969 Spin-spin interaction con- 
stants from the hyperfine structure of coupled ions Phys. Rev. 177 454-64 

[17] D'lorio IVI, Armstrong R L, and Taylor D R 1983 Longitudinal and transverse spin 
dynamics of a one-dimensional XY system studied by chlorine nuclear relaxation 
in PrCla Phys. Rev. B 27 1 664-73 

[18] Franchini F and Abanov A G 2005 Asymptotics of Toeplitz determinants and the 
emptiness formation probability for the XY spin chain J. Phys. A: Math. Gen. 38 
5069-95 

[19] Hume L and Robinson D W 1986 Return to equilibrium in the XY model J. Stat. 
Phys. 44 829^8 

[20] Lieb E, Schultz T, and IVIattis D 1961 Two soluble models of an antiferromagnetic 
chain Ann. Physics ^6 407-66 

[21] Ogata Y 2002 Diffusion of the magnetization profile in the XX model Phys. Rev. E 
66 066123-1-7 

[22] Ruelle D 2000 Natural nonequilibrium states in quantum statistical mechanics J. 
Stat. Phys. 98 57 - 75 



30 



W. H. Aschbacher 



[23] Shiroishi M, Takahashi M, and Nishiyama Y 2001 Emptiness formation probability 
for the one-dimensional isotropic XY model J. Phys. Soc. Jap. 70 3535-43 

[24] Sologubenko A V, Gianno K, Ott H R, Vietkine A, and Revcolevschi A 2001 Heat 
transport by lattice and spin excitations in the spin-chain compounds SrCu02 and 
SrsCuOs Phys. Rev. 6 64 054412-1-11 

[25] Stembridge J R 1 990 Nonintersecting paths, Pfaffians, and plane partitions Adv. 
83 96- 131 

[26] Widom H 1971 The strong Szego limit theorem for circular arcs Ind. Univ. Math. J. 
21 277-83 

[27] Yafaev D R 1992 Mathematical scattering theory: general theory (AMS: Provi- 
dence, Rl) 



